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This paper presents a conserving gapless mean-field theory for weakly interacting Bose gases. 
We first construct a mean-field Luttinger-Ward thermodynamic functional in terms of the conden- 
sate wave function "^ and the Nambu Green's function G for the quasiparticle field. Imposing its 
stationarity respect to ^ and G yields a set of equations to determine the equilibrium for general 
non-uniform systems. They have a plausible property of satisfying the Hugenholtz-Pines theorem to 
provide a gapless excitation spectrum. Also, the corresponding dynamical equations of motion obey 
various conservation laws. Thus, the present mean-field theory shares two important properties 
with the exact theory: "conserving" and "gapless." The theory is then applied to a homogeneous 
weakly interacting Bose gas with s-wave scattering length a and particle mass m to clarify its basic 
thermodynamic properties under two complementary conditions of constant density n and constant 
pressure p. The superfluid transition is predicted to be first-order because of the non-analytic na- 
ture of the order-parameter expansion near Tc inherent in Bose systems, i.e., the Landau-Ginzburg 
expansion is not possible here. The transition temperature Tc shows quite a different interaction 
dependence between the n-fixed and p-fixed cases. In the former case Tc increases from the ideal gas 
value To as Tc/Tq = 1-1-2. 33an^/^, whereas it decreases in the latter as Tc/Tq = l-3.84a(mp/27rft^)^''^. 
Temperature dependences of basic thermodynamic quantities are clarified explicitly. 



PACS numbers: 03.75.Hh, 05.30.Jp, 67.40.-w 



I. INTRODUCTION 



Mean-field theories with self-consistency equations 
have played a central role in the development of our mi- 
croscopic understanding on quantum many-particle sys- 
tems, especially on broken-symmetry phases. One of the 
most outstanding examples is undoubtedly the Bardeen- 
Cooper-Schrieffer (BCS) theory of superconductivity,^ 
which may be regarded now as the Hartree-Fock the- 
ory in the Nambu space»2, In some notable cases like the 
BCS theory, mean-field theories have brought remarkable 
quantitative descriptions of experiments. At least, each 
mean-field theory have provided a basic starting point for 
later developments, as the Weissi^ and the Stonei— theo- 
ries for localized and itinerant ferromagnets, respectively. 
Hence there is every need for a mean- field theory of every 
quantum system. However, there has been no established 
mean- field theory for Bose-Einstein condensates (BEC). 

Over the last sixty years, the field theory has 
manifested itself as one of the most powerful mi- 
croscopic approaches to quantum many-particle 
systems. It was first applied to Bose particles 
in the celebrated work of Bogoliubov,^^ followed 
by intensive theoretical investigations in later 
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However, it has brought a rather poor description of 
condensed Bose particles at finite temperatures. For 
example, we still do not have a complete agreement 
on the Bose-Einstein condensation temperature Tc of 
the homogeneous weakly interacting Bose gas under 



constant density n 



31.32.33.34.35.36.37.38.39.40.41.42.43 



not to 



mention its therm ody namic properties over < T < TcJ 
see Refs. |43 and |43 for a review on the Tc calcula- 
tions. This situation is partly due to the absence of 



a well-established renormalizcd perturbation theory 
for Bose systems corresponding to the Luttinger-Ward 
theory on Fermi liquids44 Indeed, important formal 
results on Bose systems have often been obtained 
with the bare perturbation theory, as in the case of 
the Hugenholtz-Pines theoremii^ However, the bare 
perturbation expansion itself is not very suitable for 
practical calculations due to the infrared divergences 
inherent in Bose systems. On the other hand, we still 
do not have a systematic approximation scheme of how 
to renormalize the condensate wave function and the 
quasiparticle Green's function self-consistently without 
losing the physical essentials of condensed Bose particles. 
The self-consistent Wick decomposition procedure has 
proved to be quite powerful for most quantum many- 
particle systems at finite temperatures and might be 
used also for Bose particles. The corresponding mean- 
field theory for condensed Bose systems is known now as 
the Girardeau- Arnowitt theorjii^iS^ or the Hartree-Fock- 
Bogoliubov (HFB) theoryi^ However, it predicts an en- 
ergy gap in the excitation spectrum, in contradiction to 
the Hugenholtz-Pines theorem of declaring a gapless ex- 
citation. Thus, it has become customary to introduce 
a further approximation to the HFB theory, now gener- 
ally called the "Popov" approximation;^^ of completely 
lignoring the anomalous quasiparticle pair correlation to 
recover a gapless excitation. When it is adopted to de- 
scribe dynamics, however, the same approximation does 
not satisfy various conservation laws as required. Also, 
it is not clear whether it is permissible or not to neglect 
the pair correlation completely. Recently, Proukakis et 
al.'^^ have presented an improved gapless theory with fi- 
nite pair correlation. However, it still does not satisfy 
the conservation laws. Thus, there has been no mean- 
field theory for condensed Bose systems that simultane- 



ously carries the important properties of the exact theory 
pointed out by Hohenberg and MartinsS "conserving" 
and "gapless." 

On the basis of these observations, we have recently 
formulated a new mean-field theory for BEG with the 
desired conserving gapless characteriii We here present 
a detailed description of the theory with several new re- 
sults. As first shown by Baym for normal Fermi liquids'*® 
and pointed out by Hohenberg and Martin for condensed 
Bose systemSfS^ various conservation laws are automati- 
cally satisfied in the "$ derivative approximations" where 
the irreducible self-energy E is obtained as a deriva- 
tive of some functional $ = ^[G] with respect to the 
renormalized Green's function G. This relationship be- 
tween E and G holds exactly for Fermi systems in equi- 
librium, as shown by Luttinger and Ward44 Indeed, $ 
was first introduced by Luttinger and Ward as part of 
the exact thermodynamic functional il = il[G]. It was 
later taken up by Bayni^ to present a criterion for ob- 
taining approximate dynamical equations with conserva- 
tion laws. It thus follows that dynamical equations ob- 
tained from a Luttinger- Ward functional naturally obey 
the conservation laws. We hence ask the heuristic ques- 
tion here: Gan we construct a mean-field Luttinger- Ward 
thermodynamic functional for BEG that also satisfies 
the Hugenholtz-Pines theorem? This is indeed possible 
as will be shown below. The predictions of the result- 
ing mean-field theory on the homogeneous weakly inter- 
acting Bose gas will be presented later with an expres- 
sion for ATc] this application hopefully will provide a 
deeper insight into the superfluidity and phase transition 
of still mysterious ''He^^i^ It should be noted that the 
same idea has been adopted more recently by Ivanov, 
Riek, and Knoll^"'^ (IRK) in a different context of the 
0{N) model to present an apparently different thermo- 
dynamic functional. It will be shown, however, that their 
functional yields exactly the same equilibrium thermody- 
namic properties on weakly interacting Bose gases as the 
present one. 

We finally provide three comments on terminology to 
remove possible confusions in advance. First, Gardiner— 
and Gastin and Dum^'^ have presented an alternative 
description of BEG, which was called by Gardiner— as 
"particle-number conserving Bogoliubov method." How- 
ever, his terminology is completely different from the 
present "conserving" and we call their formulations as 
"number-fixed" descriptions of BEG. Indeed, his termi- 
nology is connected more closely with the equilibrium 
thermodynamics to denote the fact that the particle num- 
ber N is chosen as an independent variable instead of 
the chemical potential /i. The two descriptions of using 
N and fi are equivalent in the thermodynamic limit ex- 
cept the fiuctuation in the particle number. It is worth 
pointing out that a "number-fixed" theory does not nec- 
essarily provide "conserving" dynamical equations in the 
present sense. Second, the term "gapless" is relevant 
here to the single-particle excitation. Indeed, the col- 
lective excitation has another story and is gapless even 



in the Girardeau-Arnowitt theory, as already shown by 
TakanoiA^ Third, it has been pointed out by Yukalo^f^ 
that the word "Popov" may not be suitable for the ap- 
proximation of setting the anomalous pair amplitude 
equal to zero in the HFB theory. Indeed, this approx- 
imation was introduced by Shohno^^ already in 1964 and 
used later by Reatto and Straleyi^-^ We can also find the 
same approximation later in the work of Baym and Grin- 
stein on the a model— We will call it as the "Shohno" 
approximation instead of "Popov" following Reatto and 
Straley.^^ 

This paper is organized as follows. Section II presents a 
mean-field free-energy functional for general non-uniform 
systems and derives a closed set of equations to determine 
the thermodynamic equilibrium. Section HI applies the 
formulation to the homogeneous weakly interacting Bose 
gas under constant density to clarify its basic thermo- 
dynamic properties. Section IV treats the same system 
under the complementary condition of constant pressure. 
Section V concludes the paper. Appendix A summarizes 
Feynman rules for the perturbation expansion with re- 
spect to the Nambu Green's function. Appendix B pro- 
vides a proof on the condensed Bose systems that the 
conservation laws are obeyed in the $ derivative approx- 
imations. In Appendix G, the connection between the 
present and the IRK theories is clarified. Appendix D 
discusses the origin of a non-Hermitian eigenvalue prob- 
lem in condensed Bose systems. Appendix E derives 
compact expressions for the equilibrium thermodynamic 
functional and the entropy. We put fcB = l throughout. 



II. MEAN-FIELD EQUATIONS 

A. Free-energy functional 

We express the field operator ?/'(r) as a sum of the 
condensate wave function ^(r) and the quasiparticle field 

(j){r) as 

V'(r) = vl/(r) + (/.(r) . (1) 

It is convenient to introduce the spinors: 



* 






$^^[<j,U], 



, #t = [^*^]_ 



(2a) 



(2b) 



Using Eq. H2a() . we define our Matsubara Green's function 
in Nambu space^ as 



OQ 

= T J2 G(r,r';z„)e- 



(3) 



with fa the third Pauh matrix and z„ = 2'KinT . The 
factor fa is usually absent in the definition of Green's 
functioniSSiii Introducing the factor brings an advantage 
that the eigenvalue problem for G~^ is equivalent to the 
Bogoliubov-de Gennes (BdG) equation for the quasipar- 
ticles, as seen below. Thus, G can be put into a diagonal 
form with respect to the eigenstates of the BdG equation. 

The bare perturbation expansion for condensed Bose 
systems may also be performed quite compactly in terms 
of the Nambu-Matsubara Green's function of Eq. ^. 
The corresponding Feynman rules can be found easily 
from those of the superconducting Fermi systems^ with 
slight modifications, which are summarized in Appendix 
A. Retaining the lowest-order diagrams in terms of the in- 
teraction and renormalizing Gq — > G, one may express the 
free-energy functional of the HFB theory in a Luttinger- 
Ward formMiS with respect to \1/ and G. 

We now present our mean-field Luttinger-Ward func- 
tional Vt = f7(\l', ^*, G) which has been obtained from 
that of the HFB theory with a slight modification. It is 
given by 

VI = [■^*{y)K^{y) rfr + - ^ Tr [ln(f3X + E-z„i) 

+ Gs]i(z„) + $. (4) 

Here K is defined by 



with 



K 



2m 



V2-|-V(r)-/i,, 



(5) 



with V(r) the external potential and /is the chemical po- 
tential, 1 is the unit matrix, and l(z„) denotes 



i(^n) 







e- 



(6) 



with 0+ an infinitesimal positive constant. The branch 
cut of the logarithm in Eq. Q is chosen along the nega- 
tive real axis, and Tr here includes an integration over 
space variables with multiplications of t^K and — z„I 
by 5{y — y') from the right implied. Finally, S denotes 
the irreducible self-energy obtained from the functional 
$ = $(*, **,G) by 



E(r,r';z„) = -- 



(5$ 



(7) 



r JG(r',r;z„)' 

With Eq. 10, Q, becomes stationary with respect to a 
variation in G satisfying Dyson's equation: 

G-^ ^Znl-hK-t. (8) 

The condensate wave function ^(r) in equilibrium is also 
determined by 5Q./5^*{y) = 0. Noting 5VL/5G = 6, we 
only need to consider the explicit ^* dependences in Q. 
to obtain 



^^^(r) = -?7(r) , 



(9) 



77(r) = 



5<^ 



(5**(r) 



(10) 



The quantity 77(r) is the so-called condensate source 
function. ^^-^ Equations Q and ifTHI constitute the <I> 
derivative approximation'^^ where the conservation laws 
are obeyed automatically by the dynamical equations. 
This is one of the main advantages in using the Luttinger- 
Ward functional as a starting point. 

Despite their importance, the conservation laws in the 
^-derivative approximation seems not to have been de- 
scribed to an enough extent in the literature for the con- 
densed Bose systems. We hence provide a detailed proof 
of them in Appendix B. 

A key quantity in Eq. Q is $. We choose it so that 
the Hugenholtz-Pines theorem is satisfied simultaneously. 
Explicitly, it is given by 



$= I dY I dY' U{y~y') 



i|*(r)nvl;(r')p 



T 



T 



^|*(r)|2Trf3G(r',r';z„)i(z„ 



- 9 E Tr f3*(r)*t(r')G(r', r; z„)i(z„) 



T 



— ^Trf3G(r,r,z„)i(z„)Trf3G(r',r';z„0i(^ 



n;n' 



T 



— ^TrG(r,r';z„)i(z„)G(r',r;z„Oi(^«') 



n,n' 



(lla) 



where U denotes the interaction potential. As may be 
realized from Appendix A where the Feynman rules in 
Nambu space are given explicitly, the five terms in the 
square brackets of Eq. Ullafl corresponds to the diagrams 
of Fig. l(a)-(e), respectively. The difference of Eq. Hlla|l 
from $HFB lies in the Fock terms, i.e., the third and 
the fifth terms. Indeed, ^hph^^^ is recovered from Eq. 
(|lla|l by replacing G and t^^ {y)^'' {y' ) in the two Fock 
terms by f3G and ^{y)^''{y'), respectively; the present 
functional was found in the reverse way. 

At this stage, it may be worth providing a comment 
on the present functional. As mentioned by Shi and 
GrifRn,M the HFB approximation yields over-counting 
of the off-diagonal self-energy diagrams in the bare per- 
turbation expansion, thereby leading to an unphysical 
energy gap in the excitation spectrum. Hence it is nec- 
essary to subtract the extra contributions from the HFB 
self-energy. In this context, it is interesting to note that 
the present functional removes some of the off-diagonal 
terms in 'I'hfb as 

$ = $HFB - /dr/"dr'Z^(r-r')[*(r)1'(r')(0^(r')0^(r)) 
-H**(r)vl/*(r')(0(r')0(r)) + |(0(r)0(r'))n . (lib) 



(a) 



(b) 



/ \ 

/ \ 

(c) 




(d) 



(e) 



FIG. 1: Diagrams contributing to $. The wavy, the sohd, and 
the broken hnes denote the interaction, the Nambu Green's 
function, and the condensate wave function, respectively. 



It satisfies the Hugenholtz-Pincs relation appropriately 
to be free from the unphysical excitation gap, as seen 
below. Thus, the subtraction mentioned above may have 
been performed appropriately. 

Based on exactly the same idea as adopted here, 
Ivanov, Riek, and Knoll^ have recently constructed an 
alternative functional $irk for the 0{N) model so as to 
satisfy the Hugenholtz-Pines theorem. Their functional 
is given in the present context by 



$ 



IRK 



$ 



HFB 



dr dr'U{r~r')\{<t){r)<t){r'))\^ . (12) 



In spite of the apparent difference between Eqs. Ijllb|) 
and l|12|l . however, $ and ^irk lead to exactly the same 
thermodynamic properties for weakly interacting Bose 
gases, as shown in Appendix C. 



B. Equilibrium solution 

Now that $ is given explicitly, we obtain the equilib- 
rium self-energy by Eq. Q. It may be written as 



I](r,r') A(r,r') 
-A*(r,r') -I]*(r,r') 



E(r,r') = 

where S and A are given by 

S(r, r') = ,5(r-r') /"z^(r-r")/9(r", r") dr' 
+ U{r-r')p{r,r') , 

A(r,r') =W(r-r')p(r,r'), 
respectively, with 

p(r,r') = *(r)M'*(r') + (0nr')'^(r)), 



(13) 



(14a) 
(14b) 
(15a) 



The equation for the condensate wave function in equi- 
librium is obtained from Eqs. ©, ll()|l . and pia|l . It can 
be written explicitly in terms of Eq. H14(l by 

K'ii{r)+ / [E(r,r')*(r')-A(r,r')**(r')]rfr' ==0. (16) 



In the homogeneous case with no external potential 
where K ^r —fx^ and VE" = ^/n^ with no the condensate 
density, Eq. (|16|) reduces to the Hugenholtz-Pines rela- 
tion /is = Ek=o— Ak=o, as desired. The expression for the 
particle number N = —dn/dfis is found similarly as 



N ■ 



/ p(r, r) dr, 



(17) 



with p given by Eq. (|15a|) . Using this equation, one 
may change independent variables of the free energy as 
n{T,V,p:)^F{T,V,N) = n{T,V,p)+pN. Equations © 
and lfT^ - lfT7|) determines the equilibrium. 



C. Bogoliubov-de Gennes equations 

Physical quantities in equilibrium may be calculated 
most conveniently in the representation where G is diag- 
onal. Noting Eq. (jH)), we introduce 



H{r,r') = T3KS{r-r') + ±{r,r') 



(18) 



The eigenvalue problem for H constitutes the BdG equa- 
tion: 



Hir,r')0^ir')dr'=E^0^ir), 



(19a) 



where (^^ denotes an eigenfunction. It has a peculiar 
feature that H is not Hermitian but satisfies ijt(r, r') = 
T^Hir' ,r)f3. Since most of the eigenvalue problems in 
quantum mechanics are Hermitian by nature, it may be 
worth tracing why we have to treat a non-Hermitian ma- 
trix in BEG. It is shown in Appendix D that this sym- 
metry necessarily results from the requirement that the 
quasiparticle field obey the Bose statistics. 

Let us enumerate basic properties of the BdG equa- 
tion, since this seems not to have been performed to an 
enough extent in the literature. We will also provide a 
compact Nambu representation to every quantity of the 
BdG equation. Taking its Hermitian conjugate and us- 
ing W{r,r') — t^H{y' ^Y)f^, Eq. (|19a|l may be written 
alternatively as 



0l,{v)T^H{vy)T:,dv^E*0l,{v'). 



(19b) 



Let us multiply Eq. H19a|) by 0h{Y)f^ from the left and 



integrate it over r. Similarly, we multiply Eq. Ijl9b|) by 
f3(^p(r') from the right and integrate it over r'. Sub- 
tracting the latter from the former, we obtain 



p(r, r') = *(r)*(r') + (0(r)0(r')) . (15b) 



{E^~EI,){0^,\t:,\0^) 



0. 



(20) 



We first put ^' = ^ in Eq. (^0)) . We then realize that E'^ 
is real as long as {0^\T^\ipfj) ^ can be satisfied by the 
eigenstate. We assume that (i) this is the case and (ii) 
eigenstates with positive eigenvalues may be normalized 
as ((/3^|f3|(/?^) = 1. We next consider the case E^^i ^ iJ^ 
in Eq. (|20(l . We then find that the two eigenstates with 
different eigenvalues are orthogonal as ((^p/|f3|(p^t) = 0. 

The Hamiltonian H has another symmetry: H{y, r') = 
— fi_ff*(r, r')fi with fi the first Pauli matrix. Let us take 
the complex conjugate of Eq. (|19a|) . multiply it by — fi 
from the left, and use the above symmetry. It is thereby 
shown that a positive eigenvalue Ei, of the BdG equation 
with the eigenfunction: 



Ui,{v) 



Uy{Y) 



is always accompanied by the negative eigenvalue 
with the eigenfunction: 



Vy{r) 



<(r) 



(21a) 
(21b) 



The sign in front of v* is introduced for convenience to 
make the coefficient Wk of the homogeneous system pos- 
itive. As mentioned above, we assume that Ui,{y) with 
El, > can be normalized as 



{u^\u^)~{v^\vi,) = l. 



(21c) 



Violation of this condition marks an instability of the 
assumed \l/(r), as is the case of a vortex- free ^(r) under 
an angular velocity f2 beyond the critical value ri^i- 

It also follows from Eq. Hl()|l that the BdG equation 
has eigenvalue whose eigenfunction is proportional to 
fsl*) with I*) given by Eq. (|2Td|) . We adopt the following 
normalized eigenfunction for this state: 



0l>%) 



1 



(^1^)1/2 



f3*(r) 



¥'o(r) 



(22a) 



It is orthogonal to \u,y) as {u^\f3\(fQ''')—0, i.e. 
{ipo\vi,). Noting Eq. (|21c|) . we may assume 

There still remains another independent state: 

¥'o(i-) 



{Ui,\^0) 



0^'\y) ^ fs0l,%) 



v'SW 



(22b) 



(22c) 



and |(^Q ) are assumed 



which does not belong to -ff, however. The functions 
^o^^'\r) (j = l,2) satisfy {^^'^\f3\4'^) = l-S,j 

The functions \uv), \v^), \(Pq 
to form a complete set. As shown by expanding an ar 
bitrary spinor with respect to the basis functions and 
calculating the corresponding expansion coefficients, this 
statement is equivalent to 



E( 



\u„)(u,y\f3 - \v^){v^\t3) 



l^0^'^)(¥?0^'V3 



\^l,'^){0. 



.'^Va 



1. 



(23) 



In summary, the BdG equation for a pair of eigenstates 
±Ei, (El, > 0) can be written compactly in a matrix form 
as 



(24) 



H{Y,Y')u,{Y')dY' ^U,{Y)E,f3., 

where iii, is defined by 



Uu{r) 



-v;{y) <(r) 



(25a) 



There exists another independent state composed of the 
condensate wave function \l/(r): 



<^o(r) 



(y9o(r) </?o(r) 



</'o(i-) 



1 



/2Nr 



:*(r), 



(25b) 

with A^o = -^(^|^) denoting the condensate number. The 
orthonormality of the basis functions reads 



f3ul{Y)T3U„,{Y)dY = S„u'l 



riV'oir)h'Po{r)dY ^ 1, 



(26a) 



(26b) 



/f3ut(r)f3(^o(r)dr = 6. (26c) 

The completeness H23|l can be written alternatively as 

Y^u,{Y)T3ut{Y')Ts + 0o{r)ri4i^')T3 ^ S{Y-Y')i , (27) 

as seen by writing down the matrix elements explicitly. 

Equations (|24|l - (|27|l are the basic properties of the 
eigenstates of the BdG equation. Noting that i^g (r) 
is absent in 7?, Eq. I|18() may be expanded as 

H(y,y') = ^7i,(r)i?,f3f3ut (r')f3 . (28) 

Green's function G= (z„l — /f)^^ becomes diagonal in 
the representation where if is. In addition, eigenvalue 
is absent in G. Hence (^(r, r'; z„) can be written as 

G(r,r';z„) = ^u,(r)G,(z„)f3ut (r')f3 , (29a) 



with 



Gv\z„ 



(z„-E,)-i 

(z„ + S,)^i 



(29b) 



Substituting Eq. if^ into Eq. © and performing sum- 
mation over Zn with t = — 0+, we obtain the expressions 
for (0t(r')^(r)) and {<^[y)<1>(y')) as 

(0t(r')<^(r)) -^ K(r)<(r')n.-K^'.(r)<(r')(l + n,)] , 

(30a) 



(30b) 
where n^= (e^"/"^— 1)^^ is the Bose distribution function, 
and we have used the (1, 2) element of Eq. H27|l to make 
Eq. (j30b|) manifestly symmetric with respect to r and r'. 
Equations l(llI)-((THJ, (EH, and ((201 form a closed set of 
self-consistent equations satisfying both the Hugenholtz- 
Pines theorem and various conservation laws. Note that 
the pair correlation ((/x/)) is adequately included in Eqs. 
(|16|l and (|24|l : neglecting this contribution yields the 
Shohno theory.^- 



D. Expressions of equilibrium 

It is shown in Appendix E that we can transform Eq. 
Q in equilibrium into 



a 



cq 



TY,Hl~e-^'^/^)-Y,E, l\v^{r)\'dr 



- /dr/rfr' 



E(r,r')p(r',r)-A*(r,r')p(r',r) 



(31) 



The entropy in equilibrium can be obtained from Sg — 
—dQ./dT by differentiating Eq. Q) in terms of the explicit 
T dependences. As shown in Appendix E, this yields 

Ss^^[{l + n^)\n{l + n^)-nA^n^] ■ (32) 
This completes the formulation of our mean-field theory. 



III. HOMOGENEOUS GAS UNDER CONSTANT 
DENSITY 

A. Equations of equilibrium 

We now apply the previous formalism to a homoge- 
neous weakly interacting Bose gas with volume V and 
particle number iV. The interaction we adopt is given by 



111 i\ 47rft2a 



(33) 



where m is the particle mass and a is the s-wave scat- 
tering length. The important dimensionless parameter of 
the system is given by 



b = ani/3 ^ 



(34) 



with n = N/V, which completely characterizes the prop- 
erties of 5<Cl. 



The corresponding condensate wave function and the 
quasiparticle eigenfunctions are the plane waves: 



*(r) = \M), 



Uk(r) 



^ikr 



V 



-Uk "k 



(35) 



(36) 



with no = Nq/V denoting the condensate density. It 
should be pointed out that the definition of ■Uk(r) above 
is slightly different from Eq. (I25all . Indeed, 'Uk(r) is com- 
posed of Uk(r) and —f^uZi^iY) with the common spatial 
dependence e*'^'", which is more convenient for the homo- 
geneous case. The coefficients Wk and Uk are chosen real; 
they clearly have the symmetry: u_k = Uk and v^k — v^. 
Accordingly, Eq. H13|l can be expanded as 



S(r,r') 



ly^gjk(r-r') 



V 



E A 
-A -E 



(37) 



where the coefficients E and A have no k dependence for 
the contact interaction. The Hartree-Fock energy S is 
obtained from Eqs. H14a() . ifTTjI . and Ip!^ as 



E = 



m 



(38) 



It merely shifts the chemical potential and does not play 
an important role in the thermodynamics with constant 
density. Substitution of Eqs. ^^ and ^7^ into Eq. l(TB|) 
yields the Hugenholtz-Pines relation: 



A*s = S - A . 



(39) 



Hence the excitation spectrum has no energy gap in our 
mean-field theory. 

The quasiparticle eigenstates are obtained by substi- 
tuting Eqs. ^ and ^ into Eqs. lO and ^ and 
diagonalizing the resulting 2x2 matrix in k space. We 
thereby arrive at the well-known expressionsi^ 



Ek^y/ek{ek + 2A), 



Uk = 



2^k 



t'k = 



Ck-^k 

2^k 



(40a) 



(40b) 



with Ek = h^k^ /2m and ^k = ek + A. Putting them back 
into Eqs. H14bp and I^Tj) with Eq. ((201, we obtain 



n^no + A ri-^.. + f A3/^ 



(41) 



A 



ATrfi^a 



m + A 



A gi/2 
E e^/T ~ 1 



de + Aey^ 



(42) 



where A = rr?^"^ l\f2'K'^'tt' is a numerical constant, and 
we have introduced an energy cutoff e^ in Eq. 1)42(1 to 
remove an ultraviolet divergence inherent in the contact 
interaction. Equations (|41|l and H42I) forms coupled self- 
consistent equations for no and A, which completely de- 
termines the thermodynamic equilibrium at a given tem- 
perature. 

Once riQ and A are obtained by Eqs. H41|l and (|42|l . 
we can calculate equilibrium thermodynamic properties 
such as the pressure Ps, the entropy S's, and the superfluid 
density ps- Pressure Ps — — ^eq/^ is transformed from 
Eq. linj into 



S A ^ [°° (lei A2 



2 2 

A 



rl/2 



iE 2E c^/T - 1 



de 



^5/2 



15\/2 
Equation l|32|) for Ss now reads 



TV 



A /-""SC + SA £3/2 



nT 



/o 3^ e 
The expression of ps is given by^ 
2A 



E/T _i 



de . 



Ps 



mn 



1 



3nT 



oo ^3/2gE/T 







(eS/T _ 1); 



de 



(43) 



(44) 



(45) 



The constant A in Eqs. 141ll - 145() under constant den- 
sity may be written alternatively in terms of the 
transition temperature of the ideal Bose gas To — 

{2nhym)[n/al)f' 



as 



A- 



2n 



V^CiDT^^' ' 



(46) 



with C(|) = 2.61 denoting the Riemann ( function. It is 
convenient to choose Tq as the unit of energy. Then Eq. 
(|1^ . for example, reads 



4e, 



1/2; 



v^cay^' 



A == 2C(|)'/3(5 



m 

n 



v^a 



A 



U/2 



E qE/t „ I 



de 



(47) 



We choose the cutoff ec so that 1 ^ Cc <C 0.3(5^^ is satisfied 
and neglect the cutoff-dependent term in the following. 
This does not affect the results qualitatively, and even 
quantitatively to the leading order in 5. 

Numerical calculations of Eqs. ((41|l - 145() can be per- 
formed easily to clarify temperature dependences of the 
basic thermodynamic quantities in the condensed phase. 
A change of variable e — 2 A sinh t is found to improve 
convergence of the integrations including e^' ^. The cor- 
responding equations for the normal state are obtained 
from Eqs. (gTJ, gSJl, and (gljl by setting A = no = and 



E — i — e + T, — Pa, which yield p^, Pn, and Sn as a func- 
tion of T. The transition temperature T^ is then deter- 
mined by the thermodynamic condition Fs{Tc,V, N) = 
Fn{Tc,V, N) appropriate under constant volume. The 
specific heat C may be calculated from the results on S 
by numerical differentiations. 



B. Properties near Tc 

We now investigate the properties near T^ (^ Tq — 1) 
in more detail. To this end, we make use of the following 
asymptotic expansions for the integrals in Eqs. (|41() - (|43|l . 
respectively: 



i 



,1/2 



V^CiDJo Ee^/T-l 



de = T3/2(i _ b^xi/2) ^ (48a) 



A 



rl/2 



V^C(i)io Ee^/T-l 



3E 



de = T^/%x^/^ , (48b) 



0?C(| 



A2\ 

2e) '^eJt 



1/2 



■ de 






(48c) 



:0.960, 6'/^C(i)/C(|) 



with X = A/r, &i ; 

1.95, and 62 = 5\/27r/6C(|) = 1-56. Each expression 
of Eq. (|48|l has been obtained as follows: (i) Expand 
the integrand in terms of A up to the order where the 
integral converges to obtain the leading terms analytic 
in X. (ii) Subtract the leading contribution from the 
original integral to express the residual term as an in- 
tegral, (iii) Estimate the higher-order integral by ex- 
panding the exponentials in the denominators up to the 
first order. For example, analytic terms are absent in 
Eq. (|48 b|) . and the integral has been estimated by ap- 
proximating c^i^ — Ik, E /T in the denominator. Those 
residual integrals generally yield terms non-analytic in 
X. This non-analyticity in the expansion near Tc is in- 
herent in Bose systems stemming from the divergence of 
the Bose distribution function at zero energy. This is 
in marked contrast to the weak-coupling theory of su- 
perconductivity where the gap equation near Tc is ana- 
lytic in x^ = (A/r)^ to make the Landau- Ginzburg ex- 
pansion possiblei^ Thus, unlike the Ginzburg-Landau 
equations for superconductorsjS the Ginzburg-Pitaevskii 
equation^ for the condensed Bose systems cannot be jus- 
tified microscopically. It is this non-analyticity which 
turns the mean-field transition into first order and also 
brings various anomalies in the thermodynamic proper- 
ties near Tc which cannot be described by the Landau 
theory of second-order transition.^'' 

Substituting Eqs. (|48a|l and (|48b|l into Eqs. ^^ and 
(I42|l . respectively, we obtain 



no 
n 



l-T^'\l-biy + b2v''). 



(49a) 



8 



no 
n 



l + h'cT^/^S. 



c5 



-Tyly 



1- 






(49b) 



with y = a;i/2 ^ (A/r)i/2, 63 = 0.559, 6^ = 1.12, and 
c = 2C(|)^/'^ = 3.79. Here we have used the resuhs 
of Eqs. H48al) and H48bl) inchiding higher-order terms 
62a; and —62a;, respectively, with 62 and 62 estimated 
numericaUyp^ Let us subtract Eq. (|49b|) from Eq. (|49a|l . 
Solving the resulting equation, we obtain y as a function 
of Tas 



viT) 



bicT^/^S 



l + (fe2 + 6'2)cTl/25 



1-- 

1/2, 



T- 



-T- 



62c,5 



(50) 



The transition temperature Tc may be estimated by set- 
ting no = in Eq. (|49b|l . which yields 



yiT, 



bicTl'^5 
l + b'cTl'^5 



0. 



(51) 



Equation (|51|l with Eq. (|50|) can be solved by expanding 
Tc in powers of 5 as 



Tr = l 



3C(i)4/3 

2.33(5+ • 



<5 + 



Stt 



3C(i)4/3 



147r 



3C(i)4/3 



-(&2 + ^)c 



(52) 



Thus, Tc initially increases linearly with 5. The numerical 
coefficient 2.33 agrees with the analytic result by Baym 
et ah?^ as well as the numerical one by Holzmann and 
KrauthiS It is worth pointing out that this expression 
makes the coefficients of 5^ and 5 in the large square 
bracket of Eq. (jSUJ) vanish at Tc. Substituting Eq. (jS^J 
into Eq. l(5T|l and noting y=(A/T)^/^, we obtain 



A(T,) 



Stt 



C(i)V3 



13.3(5-^ 



(53) 



to the leading order in 8. 

A comment is necessary on Eq. H52II at this stage. The 
expression (I52|l has been obtained as the highest temper- 
ature of the non-trivial solutions to Eq. (|49|l . However, 
the corresponding superfluid transition is first-order. It 
hence follows that, under the present conditions of fixed 
V and N , the transition temperature should have been 
determined alternatively by the requirement that the free 
energy F{T,V,N) be equal between the normal and su- 
perfluid phases. It will be shown below Eq. 1)6 5() . however, 
that Eq. 1)52(1 indeed satisfies the requirement up to the 
order 6^. 

We next consider the region below Tc and expand 
y{T) - y{Tc) in powers of l-T/T^. As noted below Eq. 
(|52|l . the leading term at Tc in the large square bracket 
of Eq. (|50|) is proportional to (5". The proportionality 



constant is found to be (fo2c)^ from the equation of order 
S^ in the expansion of Eq. H52|l . Due to this smallness of 
the constant term at T = Tc, the temperature-dependent 
contribution oc (1 — T/Tc) in the square bracket may not 
be placed outside. Indeed, we obtain 



yiT) - y{Tc) 

r0.369(5~i(l-r/Tc) 



{b,b2C^SY + ^cS{l-T/Tc, 



1/2 



1-T/Tc<^l0.56^ 



\V3C(|)'/'<5'/'(i-T/rji/2 . i-r/r,» 10.5(53' 

(54) 

Noting y=(A/T)^/^ and Eq. H53() . we may approximate 
A just below Tc as 



A(r) - A{Tc) 
(2m{l-T/Tc) 



1-T/Tc<10.5(53 



(55) 



967r (53/2 (I-T/TJ 1/2 . i_t/T,>10.5(53 

into Eq. (|49b|l . We 



Let us substitute Eqs. (|53|l and E 
then find that no{T) grows as 



'0.354(5-i(l-r/Tc) 



1-T/Tc< 10.553 



no 
n 



6tt 



.c(i)2/3 



(5^/2(i_y/2^ji/2 . i_T/Tc> 10.5(53 ' 

(56) 

Equations (|52|) . (|53|l . (|55|l . and (|55|l are the main results 
obtained from Eq. (|49|) . 

We turn our attention to other thermodynamic quan- 
tities. The chemical potential below Tc satisfies the 
Hugenholtz-Pines relation /Zs(T') = E— A(T) and directly 
reflects the singularity of A(T) near Tc- On the other 
hand, iJ-n{T) in the normal state of T>1 obeys 



C( 



3\2 



^„(r) = E-^(r3/2-i) 

47r 



(57) 



This can be shown from N = ^j^[e('^k+s-M„)/T _ j^j-i 
with the procedure of deriving the expansions of Eq. H48|l . 
Using Eqs. H52|l and (|53|) . we obtain the discontinuity 
AiJ,{Tc) = fJ,s{Tc)~Hn{Tc) as 



A^l{Tc 



An 



C(|)2/3 



^6.63(5^ 



(58) 



Next, with the procedure of deriving Eq. 1)48(1 . the en- 
tropy just below Tc is calculated from Eq. ((44)1 as 



N 



5C(I 



m) 



-T 



3/2 



-T^/^A. 



(59) 



The corresponding normal-state expression for T > 1 is 
given by 



Sn 

N 



5C(i: 
2C(i; 



2i3/2 



^C(2) /rp3/2 

87r ^ 



1) 



(60) 



Hence the discontinuity AS{Tc) = SsiTc)-Sn{Tc) is found 
by using Eqs. (|5^ and lfSS|) as 



ASiT,) 



Gtt 



c(i)2/3 



S^ = -9.94(5^ . 



(61) 



It follows from Eqs. (|55|l and H59|l that the specific heat 
Cs = T{dSs/dT) just below Tc is given to the leading 
order by 



9i 

N 



15C(f) 



9lL 

N 



4C(i) 

'4.04 



2^3/2 __ £/jn3/2 



dA 



2 dT 

l-r/rc< 10.5^3 
376^^3/2 . (62) 

i-r/rc> 10.5(53 



(l-T/T,)i/2 



It displays a divergent behavior Cg/N ^ 13.05^/^(1 — 
T/T^Y^- for 10.5(53<l-r/rc<l. 

Substituting Eq. H48c|l into Eq. ^^, ps just below T^ 
is obtained as 



Pi 



E _ A no C(i) y5/2 _ j.3/2^ 

n 2 2 n ((§) 



(63) 



The normal-state pressure of T > 1 can be calculated 
similarly to be 



p„_E , C(|)y5/2_ 



C( 



All 



(T3/2 _ 1) 



(64) 



Hence the discontinuity Ap{Tc) = ps{Tc) — Pn{Tc) is ob- 
tained as 



Ap(r,) 



47r 



C(|)2/3 



(5^ = -6.63(5^ . 



(65) 



It follows from Eqs. ^ and ^ as well as i^(r, V, N) = 
-pV+iiN that F,{Tc, V, N) = F^{T^, V, N) is satisfied for 
the expression f^ up to the order 5'^. Thus, Eq. (|^ is 
indeed a correct expression of Tc of order 5. 

We finally consider the superfluid density near Tc. 
With the procedure of deriving Eq. (|48|1 . Eq. H45|) is trans- 
formed into 



Ps 



= 1 _ 7^3/2 



4V2^rAi/2 



3C(|) 



(66) 



Substituting Eqs. (|S^ and lf5S|) into Eq. l(BB|) . we obtain 

Ap,(T,) 4^ 



mn 



3C(i)4/3 



5 = 1.16(5. 



(67) 



C. Properties at T = 

Physical properties at T = are easily found from 
Eqs. H39|l and H41|) - (|45|) by dropping terms with the Bose 
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FIG. 2: Specific heat C/N as a function of T /T^ for an^^^ = 
0.01 shown in comparison with the ideal gas result. 
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FIG. 3: Normalized condensate density no/n and superfluid 
density ps/nm as a function of T/Tc for an}''^ — 0.01. The 
ideal gas result {no/n = ps/nm) is also plotted for comparison. 



distribution function. The eading-order expressions are 
summarized as follows: 



no(0) 



= 1 ^^3/2 ^ ^ _ 1 50^3/2 ^ /gg ^ 

n 6^/71 



/i,(0) = A(0) = 2^^ = 2C(|)2/3^ = 3.795 . (68b) 



D. Numerical Results 

Figure 121 displays the specific heat CjN as a function 
of the renormalized temperature T jTc for a'n}!^ = 0.01. 
The ideal-gas result (a = 0;Tc = Tq) is also plotted for 
comparison. The curve for av}''^ = 0.01 is proportional to 
T^ at very low temperatures due to E\^(x.k for fc^O. As 
the temperature is increased, the curve starts to deviate 
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FIG. 4: Temperature dependence of the chemical potential 
H for an}'^ — 0.01 shown in comparison with the ideal gas 
result. 



upward from the ideal-gas result to show the limiting 
behavior l|^ for r<T.. 

Figure 131 plots no/n and psinm as a function of T/Tc 
for an^/'^ = 0.01 and 0. We observe that the finite interac- 
tion causes an enhancement of n^/n and Ps/nm over the 
ideal-gas value at all temperatures. The enhancement is 
larger for ps/nm than riQ/n. The quantity no/n obeys 
Eq. 1)561) near T^ and develops continuously from even 
with a finite interaction, whereas ps/nm starts from a 
positive value at Tc as Eq. (|^ . 

Figure 21 shows the chemical potential as a function of 
T/Tc for an^/^ = 0.01 and 0. With a finite interaction, 
the chemical potential displays a characteristic peak at 
Tc with a discontinuity, as shown explicitly by Eq. (|58|l . 
It eventually approaches the value of Eq. (|68bp as the 
temperature is lowered. 



E. The HFB and Shohno theories 

We now compare the above predictions of the con- 
serving gapless mean-field theory with those of the 
jjppi4,2i.29^45 ^^^ Shohno22i^ theories. Although ex- 
tensive theoretical studies have been carried out based 
on the latter theories, their predictions on the thermo- 
dynamic quantities seem not to have been clarified com- 
pletely, especially near Tc. We hence study them in detail 
with a special focus on the region T<Tc. This will also 
help us to understand common features in the mean-field 
theories of the homogenous weakly interacting Bose gas. 
The transition temperature will be determined below in 
the same way as Eq. (|52|l : the two comments below Eqs. 
(|5^ and l|H5l) also apply to the present cases. 

The equations to determine the equilibrium in the HFB 
theory are formally obtained from Eqs. H41|) and (|42|l by 
a couple of modifications: (i) a sign change for the second 
term in the square bracket of Eq. (|42|l ; (ii) the expressions 



of ^k and i?k. Indeed, they are given by 



n = no + 






A = no-T , 

m 

where F, ^k, and E'k are defined by 



■E 



k 



(69a) 



(69b) 



(70) 



Ck = ek+A+2r, Sk= V(ek + 2r)(ek + 2A + 2F), (71) 

with e)s_^h^k^ /2ra. The Shohno theory is obtained from 
Eqs. ^ and (ETJ by setting F = 0. 

Near T^Tq — 1 where A and F are much smaller than 
To, Eq. 1)69)1 is approximated to the leading order by 



no 
n 



= 1 - T^'^- [l - bi (x/F+A + Vf )] , 



c6 



no 
n 



6i(Vr+A-Vf) 



(72a) 



(72b) 



with bi = \/2^/C(i) and c = 2C,{%fl'^. These equa- 
tions correspond to Eq. ()49)l and have been obtained 
with the procedure of deriving the expansions of Eq. 
()48)l . Note that Eq. ()72)l is also non-analytic in F and 
F-t-A. Using A= c(5(no/n — F) from Eq. ()69b)l . we next 
rewrite Eq. I)72b)) in terms of F and no/n. We also in- 
troduce the dimensionless variables r, s, and t through 
\/ riQ / n = bi^/cS r , vr = &iC(5s, and T=l-\-b\c5t. Retain- 
ing terms of order 5, Eq. 1)72)1 is transformed into 



2 3 

Equation (|75|) yields 



-l + \/TT4f 



(73) 



f{r,t) 



2r^ + {2+it)r^ - {2+it)r + -t{2+3t) = . 



(74) 
It has a couple of solutions ri , r2 > for t w but no 
real solution for a large-enough t. Thus, the transition 
temperature corresponds to the point where / = and 
df /dr — 0, which yields the critical values of i = 11/24 
and r — 3/4. Thus, the transition temperature of the 
HFB theory for (5->0 is given by T^^^ = 1+AT^^^ with 



AT 



HFB 



IItt 



6C(i)*/^ 



(5 = 1.60(5. 



(75a) 



The transition is clearly first order. The discontinuities 
of various thermodynamic quantities at T^^^ can be ob- 
tained with the procedure of Sec. lIIIBl as 



,HFB 



(^c) 



97r 



4C(|)4/3 



1.96(5, 



(75b) 
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, HFB 



[Tc 



Stt 



2C(i)2/3 



6^ = 4.14(5^ 



(75c) 



p»^«(r,) 139^ 



mn 60C(|)'*/3 



(5 = 2.02(5. (75d) 



Note that even no/n is discontinuous at Tc in the HFB 
theory. 

The same analysis can be performed for the Shohno 
theory. We thereby obtain 



AT; 



'Shohno 



27r 



3C(i)^/3 



j^bhohno^jnj _ TT 



C(|)4/3 



(5=^0.582(5, (76a) 



5 = 0.873(5, (76b) 



A'"'°""°(Tc) 



27r 



C(|)2/3 



5^ = 3.31(5^ , (76c) 



^Shohno 



{Tc) 



3C(|j^' 



1.46(5. 



(76d) 



Thus, even the Shohno theory predicts an enhancement 
of Tc over the ideal-gas value Tq. 

As for the overall temperature dependences of various 
thermodynamic quantities, the HFB and Shohno theories 
both yield results qualitatively similar to those of the 
previous mean-field theory except the discontinuity of tiq 
at Tc- For example, the specific heat of the HFB and 
Shohno theories also display the divergent behavior of 
(rc-r)-i/2 j^s^ ^qIq^ rj.^^ o^g ^Q a gj^itj; Y, the HFB 

theory generally predicts larger deviations from the ideal- 
gas results than the Shohno theory. 



IV. HOMOGENEOUS GAS UNDER CONSTANT 
DENSITY 

We finally study the homogeneous weakly interacting 
Bose gas under the complementary condition of constant 
pressure. Despite its fundamental importance as a sub- 
ject of quantum statistical mechanics, the system seems 
to have been investigated only by Reatto and StralejsSS 
based on the HFB and Shohno theories iS2ii^ Moreover, 
their main interest was on superfluid "^He so that the 
properties in the weak-coupling region remain essentially 
unexplored. 

Let us introduce a dimensionless parameter 5p by 






(77) 



which completely characterizes the weak-coupling region 
of (5p <C 1. The thermodynamic equilibrium of the con- 
densed phase is determined by Eqs. I@lj-(@3l- They yield 
the density n, the condensate density riQ, and the pair po- 
tential A as a function of p and T . The prefactor A in 
Eqs. (|41|l - l|43() may be rewritten conveniently as 



A^ 



2nr. 



2p 



v^a 



1)T 



3/2 



v^ai)To' 



(78) 



where ric and Tq are the critical density and the transition 
temperature of the ideal Bose gas, respectively. They 
are given in terms of p as Uc = [mp/27r/i^^(|)]'^/^(^(|) 
and To = (27r?iVTO)3/5[p/(^(|)]2/5. The equations for the 
normal state are obtained from Eqs. I|41(l . H43(l . and H44|l 
by setting A = ng = and E = £^ = e + Sn — ^n; they 
yield n, ii^, and ^n as a function of p and T. The transi- 
tion temperature Tc is then determined by the thermody- 
namic condition ^s{Tc,p) = finiTcp) appropriate under 
constant pressure. The specific heat per particle Cp may 
be calculated from the results on S by numerical differ- 
entiations. In the following, we will adopt the units in 
which nc=TQ = l and neglect the cutoff-dependent term 
in Eq. (|42|l as was done in the case of constant density. 

To investigate the properties near Tc (^1) in more 
detail, we make use of the asymptotic expansions of Eq. 
(|1S|| . Equations ^^ and ^^ with Eq. |||SJ) enable us to 
express uq and n in terms of T and A as no = A/c'dp — 
&1TAI/2 and n = T^^^ + A/c'Sp - 2biTA^/^ with c' = 
2C(|)/C(|)^^^ = 4.93. Substituting them into it, Eq. lO 
becomes a nonlinear equation for A as 



b''] TA3/2 



b'l ,2 _ (JhK 

2c'dp V 2 

+ (6'/r3/2+4fe2^;'c'r2(5p) A - 46i6'/c'T^/2(5pAi/2 

^^5/2 _^ 1,'ic'T^Sp -1 = 0. (79) 

Thus, the equation is non-analytic in A and completely 
different from the one in the Landau theory of second- 
order transitioniffi It is convenient to rewrite T in Eq. 
(El as T = 1 - f 6" c'(5p + tSl = l- 3MSp 
is a parameter appropriate to describe the region T 
Then Eq. H79|l can be solved analytically by expanding A 

as A = A2it)6l+ A3{t)S^-\ . The expression for A2{t) 

is obtained as 



tSp, where t 
' ~ 1. 



A2(i) 



26ic' + J(26ic')2 



9b'{c'' 
10 



26'/ 



(80a) 



It hence follows that a superfluid solution exists for t < 
U = {2b'{/5)[{2bic')^ + 9b'{c''^/10] = 103. The quantity 
A3 (t) may be expressed in terms of t and A2 (t) . A similar 
consideration for the normal state leads to a nonlinear 
equationfor S— /i. Expanding!]—^ = T,2{t)Sp+'S3{t)S'p + 
• • • , S2 (0 is found to satisfy 



VS^=-x/26ic' + J2(6ic')2 



9b'{c'' 
10 



2b'l 



t. (80b) 
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FIG. 5: Specific heat Cp as a function of T/T^ for Jp = 0.001 
and 0.005. 



FIG. 7; Chemical potential as a function of T/Tc for Sp ■ 
0.001 and 0.005. 
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FIG. 6: Normalized condensate density no/n as a function 
of T/Tc for 5p = 0.001 and 0.005. 



FIG. 8: Normalized particle density n/ric as a function of 
T/Tc for 5p = 0.001 and 0.005. 



It has a real solution E2 for t>t„ = 18(6'/c')V50 = 33.1. 

From ts>tn, we realize that the superfluid transition 
is first order with a metastable region extending over 
tn ^ t ^ts- The transition temperature T^ can also be 
expressed as 



Tr 



1 - 3.84,5p + tj^ . 



(81) 



Thus, Tc becomes smaller than To = 1 to the leading order 
in Sp. This behavior is opposite to the case of constant 
volume where Tc increases as a function of the dimen- 
sionless parameter an^/'^.^^ Indeed, the effect of the weak 
repulsive interaction is completely different between the 
volume-fixed and pressure-fixed cases. In the former case, 
a finite interaction suppresses the density fluctuation so 
that it works favorably for the phase coherence over the 
system, thereby leading to an enhancement of Tc- In the 
latter case, on the other hand, the interaction lowers the 
particle density to decrease Tc- The value of the pref- 
actor 3.84 is almost certainly correct, in contrast to the 
volume-fixed case where it seems still controversial42*^ 
We turn our attention to tc of Eq. (|81|l which includes 
all the higher-order contributions. It is found that tc 



has a rather large Sp dependence, due partly to the large 
values of the numerical constants appearing in the equa- 
tions, e.g., is = 103. It decreasesfrom tc((5p— >0)=ts = 103 
down to 98.6, 96.9, and 95.7 at Sp = 5.0 x lO"'^, 1.0 x IQ-^, 
and 1.5 x 10~^, respectively. The contributions from A3 
and S3 are essential to produce the results, which agree 
excellently with the exact numerical results of using Eqs. 
(|41|l - (|43|) for T ~ 1. Discontinuities in the thermody- 
namic quantities at Tc are found to be all of the order of 
Sp in accordance with the tiny metastable region. 

Let us move onto the overall temperature dependences 
of various thermodynamic quantities. Figure 5 plots spe- 
cific heat Cp over 0.6Tc < T < 1.2Tc for Sp = 0.001 
and 0.005. A remarkable enhancement of Cp near Tc is 
clearly seen. Indeed, the curves show a divergent behav- 
ior (x\T—Tc\~^^^ on both sides near Tc- It is completely 
different from the discontinuous behavior predicted by 
the Landau theory of second-order transition and caused 
by the non-analytic nature of the order-parameter ex- 
pansion near Tc- The sharp heat absorption/emission 
near Tc may be realized in terms of the proximity to 
the ideal Bose gas under constant pressure which sud- 
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denly falls into the ground state of volume zero below 
Tq. The weak repulsive interaction weakens the singu- 
larity of the ideal Bose gas at T = Tq over a finite range 
close to Tc- This sharp absorption/emission of heat near 
Tc in the present system is in qualitative agreement with 
the observation on ^He4^ Figure 6 shows the normal- 
ized condensate density no/n as a function of temper- 
ature, where a steep change just below Tc can also be 
seen clearly. Although not plotted here, the normal- 
ized superfluid density Ps{T)/mn{T) is almost indistin- 
guishable in the present scale from no{T)/n{T) for both 
i5p = 0.001 and 0.005. Figure 7 displays temperature de- 
pendence of the chemical potential fj,. The curves flatten 
out rapidly below Tc reflecting proximity to the ground 
state. Finally, Fig. 8 plots the particle density n as a 
function of temperature for 6p = 0.001 and 0.005. Since 
the superfluid component contributes less to the pres- 
sure than the normal component in the weak-coupling 
region, a constant pressure can only be sustained at 
lower temperatures by increasing the particle density. 
Thus, n{T) increases monotonically as T is lowered from 
Tc- The leading-order expression of n(0) is given by 

n(0)/nc = [C(|)'/VC(i)]<^p"'^' = 0.457<5p-'/^ implying a 
smaller enhancement of n{0) for a stronger interaction. 
In this respect, it is worth pointing out that n{T) of 
superfluid ''He displays a maximum at Tc^ Whether 
n{0)/n(Tc) may be reduced beyond 1 or not in the strong- 
coupling region remains a challenging problem to be clar- 
ified theoretically. 



V. SUMMARY 

We have constructed a mean-field theory for Bosc- 
Einstcin condensates. The basic equations to determine 
the equilibrium are given by Eqs. (|14|I - H18|) . (|24|l . and 
(|30|l . They can be applicable to general non- uniform sys- 
tems such as the trapped atomic gases with or without 
vortices. The present theory has a manifest advantage 
over the HFB and Shohno theories that it simultane- 
ously satisfies the Hugenholtz-Pines theorem and vari- 
ous conservation laws. It will be useful in clarifying both 
equilibrium and dynamical properties of condensed Bose 
gases over the whole temperature range. The conserving 
gapless theory is expected to provide at least a qualita- 
tive description of weakly interacting BEC. Especially, it 
will shed a new light on the temperature dependences of 
the collective modes in BEC. It may also form a starting 
point for a unified description of BEC over the whole in- 
teraction strengths from the trapped weakly interacting 
atomic gases to superfluid *He. 

The theory is then applied to a homogeneous weakly 
interacting Bose gas with s-wave scattering length a un- 
der constant density n in Sec. Ill and under constant 
pressure p in Sec. IV. The order-parameter expansion 
near Tc is found not of the Landau-Ginzburg type but 
non-analytic by nature due to the divergence of the Bose 
distribution function at zero energy, as shown explicitly 



in Sec. IIIB. This non-analyticity makes the superfluid 
transition first-order and also brings unique temperature 
dependences of various thermodynamic quantities near 
Tc which are quite different from those of the Landau 
theory of second-order transition. The transition tem- 
perature Tc under constant n is found to increase from 
the ideal Bose gas value Tq as Tc/Tq = l + 2.33an^/^To, 
whereas it decreases under constant pressure as Tc/Tq — 
l—3.84:a{mp/2'Kh'^y^^ We have also clarified overall tem- 
perature dependences of basic thermodynamic quantities 
in Sees. HID and IV. These predictions are then com- 
pared with tho se of the BFBidiSliSSJk and ShohnoSSiSa 
theories in Sec. IIII El Those theories also predict a first- 
order transition and an enhancement of Tc- 

We finally comment on the first-order superfluid tran- 
sition of the present mean-fleld theory. As mentioned 
above, it is a feature inherent in the mean-field theo- 
ries of BEC stemming from the non-analyticity of the 
Bose distribution funtion at zero energy. However, it is in 
contradiction to the observation of superfluid ^He where 
the transition is continuous. It has been pointed out 
by Reatto and Straley^S that the mean-field transition 
would be continuous if the change of the single-particle 
dispersion at Tc were from fc*^ to fc' " with e < f {k: wave 
number) . The behavior fc' with e < § above Tc may only 
be possible by including correlations and/or fluctuations. 
Thus, we encounter the question of whether the correla- 
tions and/or fluctuations are essential even in the limit 
a -^^ to transform the nature of the superfluid transi- 
tion. If not so, the transition is expected to change its 
character from flrst-order into continuous at an interme- 
diate value of a. Thus, the predictions of the mean-field 
theory raise the questions on: (i) the relative importance 
of the single-particle and collective excitations in BEC as 
a function of the interaction strength; (ii) the role of the 
single-particle excitations in superfluid '*He above and be- 
low Tc- It is worth reminding that the transition of the 
ideal Bose gas under constant pressure is strongly first- 
order. Hence it may not be so unreasonable to expect 
that the first-order transition persists into the region of 
a finite interaction. 

Thus far, microscopic theoretical studies on condensed 
Bose systems seem to have been carried out sepa- 
rately in two limiting cases with completely different 
approaches: (a) the weak-coupling regime where cer- 
tain model Hamiltonians have been investigated based 
on the field theory^^S*^ (b) the strong coupling regime 
where interests have been focused on describing su- 
perfluid *He quantitatively based on some variational 
wave functions r2iSiS4 lattice-gas models which can be 
mapped onto spin systemsjS and quantum Monte-Carlo 
calculations.^^ We hence need systematic investigations 
based on a model Hamiltonian from weak- to strong- 
coupling regimes, which may provide us more profound 
understanding on condensed Bose systems and shed a 
new light on the nature of the still mysterious lambda 
transition of superfluid "'He.^ 
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APPENDIX A: FEYNMAN RULES IN NAMBU 
SPACE 

The bare perturbation expansion of the condensed 
Bose systems is usually carried out with respect to the 
condensate wave function ^(r) and the non- interacting 
Green's function: 

G(°)(riri,r2T2) = -(r,(/.(riTi)0^(r2r2))o. (Al) 



Ct 



rule: Contract ipi only with (j)j in ^4'j4'j ^-nd multiply 
the contribution by 2. This yields (j)UTrN(l>j(j)J)Q, where 
N is effective only when i and j belong to the same inter- 
action line according to the definition in Eq. ljA3p . Note 

{Tr$j4>})o = -TzG^^^ as realized from Eq. ©. The rule 
allows us to perform a perturbation expansion in Nambu 
space by using only the Feynman diagrams for G- . The 
following two points should be supplemented: (a) A fac- 
tor 2 remains uncancelled for each closed particle line, 
due to the fact that we have to fix the starting pair 5 </'/0i 
for each closed line, (b) The contraction of 0^ with the 
last (j), in the closed particle line may be transformed as 



^{TrN$lAk^$^)o = ^Trife, (T, A^0,,?,t)o , (A5) 



The anomalous (i.e., off-diagonal) Green's function is not 
necessary before the renormalization, and we have well- 
established Feynman rules for this purpose. We will show 
that the same rules can be used with a slight modification 
to the perturbation expansion in terms of the Nambu 
matrix of Eq. Q. This extension will be performed in 
exactly the same way as the extension from the normal 
to the superconducting Fermi systems. '"^^ 

Let us specifically consider the following term in the 
interaction: 



i?(3) = [/^^,(^*0t,0^,^ 



i + 4>\(l)\,<i3v^i 



(A2) 



(3) = UivN(^^14>^4>\,^v + 4>\^i4>\'4>i) ■ (A3a) 



where Z^ii'=Z^(ri—r'^), (j)i = (j){vi), etc., and summations 
over the repeated indices are implied. Using the normal- 
ordering operator TV, we first rewrite Eq. (jA2|l as 



H 



With the operator A^ in front, Eq. (|A2p is further trans- 
formed into 

(A3b) 
where 4> and ^ are the spinors of Eq. |5J. Comparing 
Eqs. l)A3a|l and (|A3b(l . we observe the following corre- 
spondence: 



where A^^ is a Nambu matrix composed of the contrac- 
tions connecting k with i. (c) As for the open line start- 
ing from ^(pl'^i and ending at ^^^0^, the two factors 
2 remain uncancelled. These rules are valid even after 
including the interactions other than Eq. IJA3|I . 

In summary, we can carry out a perturbation expan- 
sion in Nambu space by modifying the Feynman rules for 
G(o), *, and ^* as follows: (i) GJf ^fgClf, ^^^\^^, 
and ^* — > \'^i, (ii) ^Tr for every closed particle line. 
Note that the operator N in {TrN(l)i(f>l)o transforms into 
the matrix of Eq. © in the frequency space. 

Let us write down the contributions of the diagrams 
in Fig. n based on the above rules with G^°^ -^ G. This 
yields <1>hfb- We then modify it slightly as described 
below Eq. (|lla|l . We thereby obtain (f> of Eq. (|lla|l . 



APPENDIX B: CONSERVATION LAWS 

It is shown that the $ derivative approximation de- 
fined by Eqs. Q and HlOf) automatically satisfies various 
conservation laws. 



Definitions 



?^ 



*i ^ * 



t 
1 ' 



4^1 



61 , etc. (A4) 



Let us perform a perturbation expansion in the imaginary 
time domain with respect to (|A3b|l . We hence include t^ 
in the index i, i.e., i = riTi and z' = r^rj. Consider an nth- 
order term in the expansion and take a contraction of (/)/ 



with i 



:t. 



,- . We then find that $] and 



con- 



.j V'j- "^ "^i^ii ii^i^ """" Yj ""^ Yj ^" 2l'j'f'3 

tribute equivalently to the contraction, as shown easily 
by writing it down with respect to 4>i, (j)], 4>j , and (/)] ; the 
statement is valid even if we formally retain the anoma- 
lous averages. Because of this fact, we can introduce the 



Since we work in the real-time domain, we need mod- 
ifications of various definitions. The Green's function is 
now defined by 



with l = riti. It obeys the equation of motion: 
_9G(1,1') 



(Bl) 



ih- 



dti 



f3XiG(l, 1') - 1<5(1, l') + S(l, 2)G(2, 1') , 

(B2a) 



15 



or equivalently, 
_dG{l,l') 



-ih- 



dtv 



ifi,G(l, l')f3 = W(l, l') + G(l, 2)S(2, 1') 



(B2b) 

Here Ki is given by Eq. |Hll, and integration over the 
barred index 2 is implied. Next, the condensate wave 
function satisfies 



oil 



or equivalently, 
d 



_j?i *t(i) _ Ki¥{1)t3 = r7t(l)f3 
oil 



Here ^ is given in Eq. jSbJ), and ^(1) denotes 

,7(1) 



^7(1) 



^*(1) 



(B2c) 



(B2d) 



(B3) 



The key quantities in Eq. IJB2|) are S and rj. In the $- 
derivative approximation^'^ adopted here, they are deter- 
mined from a functional $ = $((?, ^, ^*) by 

(5$ 



2(1,1') = -2i 



5G(1',1) 



77(1) - r^ 



5$ 



<5**(1) 



(B4a) 



(B4b) 



Note that this $ is different in definition from those of 
Bayn>^* and Hohenberg and MartinS^ by factor i, i.e., 
4> — i<&Baym = ^'I'hm- The present definition has an ad- 
vantage that $ becomes real within the mean-field ap- 
proximation. A comment is also necessary on Eq. (|Bip . 
When t2 — ti, the operator (p^ in G(l, 2) should be placed 
to the left of 4> by definition. This procedure is in accor- 
dance with the treatment on the equal-time operators in 
the calculation of $. Put it another way, the equal time 
implies ^2 =ti+ for the (11) component of G and ^2 = ii- 
for the (22) component. Thus, ^2 is not exactly equal to 
ti for the diagonal elements of G. 



2. Identites 

With Eq. (|B4p . we can prove several identities by ex- 
actly following the argument of Baym4& 

Suppose we perform the gauge transformation: 



'g(1,2) -^ e^x(i)^3G'(i^2)e^*x(2)7 
*(1) ^e*x(i)*(l) 



(B5) 



where x(l) is an arbitrary real function. The correspond- 
ing variations in G and '5 are given to the leading order 
by 

r<5G(l, 2) = *[x(l)f3G(l, 2) - G(l, 2)x(2)f3] 
\5*(1)=*X(1)*(1) ' ^ ' 



respectively. However, $ is clearly invariant under Eq. 
(|B5|) . i.e., (5$ = 0. This yields 

yTrS(T, 2)5G(2, 1) + ,5**(1)77(1) + S^{l)rj*{l) = . 

(B7) 
Substituting Eq. HB6|) into Eq. HB7|) and recalling x(l) is 
arbitrary, we obtain 

iTrf3[S(l,2)G(2,l)-G(l,2)S(2,l)]-^^t(i)^3^(l)^0, 

(B8) 
with *+ = (***). 

We next consider the Galilean transformation: 

fG(l,2)->G(ri-|-R(ti),ti;r2-fR(t2),t2) 
|*(l)^*(ri+R(ii),ti) 



(B9) 



The corresponding leading-order variations in G and ^ 
are given by 

f5G(l,2) = [R(fi).Vi+R(f2)-V2]G(l,2) 
\<5vl/(l) = R(fi).ViVl/(l) 

However, $ remains invariant since it is a quantity ob- 
tained by integrations over all the space-time variables. 
We hence conclude that Eq. IjBZp holds also in this case. 
Substituting Eq. IjBlOp into Eq. (|B7|I and recalling R(i) 
is arbitrary, we obtain 



fs{l)dri =0. 



(Bll) 



with 



S(l) = -(Vi- ViO{*ftTr[G(l, 2)E(2, 1') - S(l, 2) 



xG(2,l')]-2*t(i')^-(l)}|^,^^, 



(B12) 



Here terms with ViS(l,2) and Vi77'(l) have been ob- 
tained through integrations by parts. 

We finally consider the "rubber-clock" transformation 
oi Baym-^ t^e{t) = t+ fit) with /(±oo) = 0. It yields 



-(^^^)-£r-(-^-^^^)(sr 



*(i) 



d9i_ 
'dh 



1/4 



(B13) 



*(ri0i) 



The factor {dO/dt)^/'^ compensates the Jacobian for t—tO, 
thereby making the integration of H^t (t) over t invariant 
in form through the change of variables. The correspond- 
ing leading-order variations in G and VP are given by 



5G(1,2) 



^ + < + /^J^)^(^'^) 



5*(i)=(4+/i|^)*(i) 



(B14) 
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with fj = f{tj). However, <I> remains invariant. Hence 
Eq. (|_B7|I holds also in this case. Substituting Eq. (|1J14() 
into Eq. I|B7|I and recalling f{t) is arbitrary, we obtain 
the third identity: 

^^^ = \ l^ri A{zOTr[G(l, 2)S(2, l') + S(l', 2) 
xG(2, 1)] +*t(l),7(l') + ?7^(l')*(l)}i,=i , (B15) 
where (i?int)i is defined by 

(iJi„t)i = ^ /"dri{iftTr[G(l,2)E(2,l) + E(l,2)G(2,l)] 
+^t(i),y(i)+^t(i)^(l)}. (B16) 



Equation ljB16p indeed corresponds to the interaction 
energy of the system. This can be checked as follows: 
(i) write down the equations of motion for and ^ in 
terms of the field operators and the condensate wave 
function; (ii) construct the interaction energy by adding 
those equations appropriately; and (iii) transform the 
one-particle energy in the equation by using Eq. (|B2|I . 

Equations (|B8(I . (|B11(1 . and l)B15|l are the basic iden- 
tities to be used below. 



3. Conservation laws 

Let us operate ^Trfs , — iXrfa, and — ■|^^(1) to Eqs. 
(|B2ap . (jB2bp . and iJB2c|) . respectively. Adding the result- 
ing three equations, setting 1' = 1, and noting Eq. (|B8p . 
we obtain the particle conservation law as 



dnjl) 



Vij(l)=0, 



where n(l) and j(l) are defined by 
n(l)=p(l,l). 



(B17) 



(B18) 



Here j(l), n(l), and S(l) are given by Eqs. ||BT9)| . ||BT8)| . 
and 1B12II . respectively, and n(l) is a tensor defined by 



n„;3(i) 



Am ' 



;(Vi„-Vr„)(Vi;3-Vr;5)p(l, l')|i,^i ■ 

(B22) 

Integrating Eq. IJB21|I over ri and noting Eq. l|Blip . we 
obtain the conservation law for the total momentum as 



_d_ 



mdvi 



n{l)ViV{l)dri 



(B23) 



We finally consider the energy conservation law. The 
time derivative of the interaction energy is given by Eq. 
(jBlSp . The terms in the curly bracket of Eq. IjBlSp can 
be transformed into an expression of using the left-hand 
sides of Eq. (IB2p . We then find that the terms with the 
operator ifi-gfy cancels out upon applying -^ in front of 
the curly bracket. We thereby obtain 



d{Hn 



drii^i— p(riii,r'iti)|ri=ri , (B24) 



dh 



with p(l, 1') defined by Eq. l|B20)l . Wh en K has no 
time dependence as given by Eq. lO, Eq. I|B24|I directly 
leads to {K + Hint) = constant. Following Kadanoff and 
Bayruji^ we may alternatively consider y as a small time- 
dependent external field and trace the time variation of 
the internal ener gy {K^^^+Hi^t), where K'^^^ = -^V^. 



Using Eq. ljB17p and performing integration by parts 
with respect to ri, we obtain 



d{K(°')+Hn 
dti 



- I i{l)-ViV{l)dri. (B25) 



Thus, we have seen that the conservation laws are au- 
tomatically obeyed as Eq. (|B17|) . (|B23I) . and HB25|I in the 
<& derivative approximation of Eq. IIB4I) . 



APPENDIX C: THE IRK THEORY 



j(l) 



—ih 
'2m 



(Vi-ViOp(1,1')|i,.i 



with 



p(l,l')^(0Ul')0(l)>+*(l)**(l'). 



(B19) 



(B20) 



This p{l, 1') is the one-body density matrix appearing in 
the evaluation of every one-particle operator. 

To obtain the momentum conservation law, we operate 
-£Tr(Vi-ViO, i|Tr(Vi-ViO, and -^^\l)f,V, to 
Eqs. (|B2a|) . (|B2b|l . and (|B2c|) . respectively. We also mul- 
tiply Eq. (|B2d|l from the right by 2^f3Vi*(l). Adding 
the four equations and setting l' = l, we obtain 



^J(1) + ^V,F(1) 
oil fn 



-Vi-n(l) = — S(l). (B21) 



It will be shown that the thermodynamic functional 
presented by Ivanov, Rick, and Knoll^^ (IRK) leads to ex- 
actly the same thermodynamic properties for the weakly 
interacting Bose gases as the present theory. 

The IRK self-energy s(ii^k) jg calculated from Eq. IpH) 
by Eq. {71). The result can be written as 



f](IRK)(^ 



r,r = 



I](r,r') „A(iRK)(r,r') 



(CI) 



where S is given by Eq. (|14a|l . and A^^^^^ denotes 

A(^^^\r,r')=U{r-r')p'^^^^\r,r'), (C2) 

with 

^(IRK) (j.^ J.,) ^ *(r)*(r') - (0(r)0(r'))(™^) . (C3) 
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The minus sign in front of A^^^^^ in Eq. IJC1|I has been where 

introduced to make the coefficient A^j^ ' for the homoge- 
neous system positive at low temperatures, in accordance 
with the definition of A in Eq. (|14b|l . Next, the equation 
for ^(r) is obtained from the thermodynamic functional -j- _ 

^IRK ^jy dQ}^^/d-^{r)^Q. It is given by Eq. ^ with ^ = 

the replacement: A^A^^"^'. 

In a way parallel to Eqs. IjlSI) and H24() . let us introduce 
the Hamiltonian: 



^(IRK) (j.^ J./) ^ f3i^j(r _ r') + i](iRK) (j.^ J./) ^ ^C4) 

and construct the eigenvalue problem as 



/'^(IRK)(j. j.,)^(IRK) 



(r')dr' 



-(IRK) 



where v/u \s defined by 



-(IRK) 



(r) 



vI{y) <(r) 



{r)E,T^ , (C5) 



(C6) 



The sign of v*{y) is chosen so that the coefficient Uk 
of the homogeneous system is positive, in accordance 
with the definition of w*(r) in Eq. 125a|) . The corre- 
sponding Green's function G^^^ is obtained from Eq. 
(|29|l with the replacement u^ ^ ul, '. Noting Eq. 
(j30bp and the difference between Eqs. (|25a|l and ljC6p . we 
now conclude that ((/)(r)(/)(r'))(^^^) in Eq. IjCSp satisfies 
-(0(r)0(r'))(^^'^) = (0(r)0(r')), which implies A(iRK) ^ 
A. Thus, the equations to determine the equilibrium are 
exactly the same between the present theory and the IRK 
theory. Finally, one can show that the equilibrium ther- 
modynamic potential of the IRK theory is also given by 
Eq. |j23l. Thus, the two functionals Hlla|) and ((T^ yield 
exactly the same thermodynamic properties on BEG. 



APPENDIX D: NON-HERMITIAN EIGENVALUE 
PROBLEM 

In this Appendix we study possible origins of the non- 
Hermitian eigenvalue problem in BEG. It is shown that 
we have to treat the non-Hermitian matrix of Eq. H24|l 
necessarily if we impose the two requirements: (i) the 
quasiparticle field obey the Bose statistics; (ii) the quasi- 
particle eigenstates be obtained by diagonalizing a ma- 
trix. 

Let us introduce the quasiparticle field ^^ through 



E 






rv 






-ll 



(Dl) 



where v distinguishes quasiparticle eigenstates, and 0r = 
0(r), Uri, — u^{r), etc. Equation (|D1I) may be written 
symbolically in a compact form as 



U7, 



(D2) 



- <^r " 








" 7i^ " 


0r' 








liy' 


-0t 


, U = 


U —V 

-v*u*_ 


, 7 = 


-ll 













(D3) 



with 



{u)yy 



ilDr 



We now impose the condition that j^, as well as 
the Bose statistics. This implies 

where t^ is an infinite matrix defined by 

1 



T3 



-1 



(D4) 
obey 

(D5) 
(D6) 



with 1 and denoting the unit and the zero matrices, re- 
spectively. We also require that u be composed of eigen- 
states of a matrix H as 



Hu — uE , 



(D7) 



where i? is a real diagonal matrix. Multiplying Eq. IjDTI) 
by u^^ — fsM^fa from the left, we obtain 



u)hHu^E. 



(D8a) 



Tsu'Ta 

We next take the Hermitian conjugate of Eq. (|D8a|) and 
multiply the resulting equation by fa from both sides; 
this yields 

f3u''H^T3U = E. (D8b) 

Subtracting Eq. ljD8b|) from Eq. ljD8a|) and multiplying 
the resulting equation by T3U and t^u^ from the left and 
the right, respectively, we obtain 



ijt ^ fgijf. 



(D9) 



Thus, H is necessarily non-Hermitian with the property 
(lD9ll. 



APPENDIX E: DERIVING EQS. ^^ AND STUi 

We here derive expressions of the equilibrium thermo- 
dynamic potential H31I) and the entropy H32|) . 

Using Eqs. H18|l and (|28|l . the logarithmic term in Eq. 
Q can be transformed as 

L = ^ ^ Tr i(z„) \n{T3K+±~zJ) 

n 

J Jc ^ 

+ ln(z-^,)] - |w,(r)|2[ln(£;, + z)+ln(-z-£;,)]} , 



(El) 
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with ri(z) = l/(e^/"^— 1). Here we have used the standard 
procedure to change the summation over n into a contour 
integrals^ with the replacements e^"°+ -^ n(2;)e^°+ and 
e-2nO+ _^ _,^(_2;)e~^°+ in the matrix i{zn) of Eq. ^. 
Contour C consists of two parallel lines run ning along 
the imaginary axis as shown in Fig. 25.4 of Ref.|63, which 
can be deformed into contour C" in the same figure along 
the real axis. Performing integration by parts in terms 
of z and collecting residues within C", we obtain 



L = Y, T\n{l-e-''"/^)-E^ f\v,ir)fdr 



(E2) 



Next, we rewrite $ of Eq. Hlla() by using Eqs. l(T ^ - l(T5|) 



as 



$ = -Tr I dr /dr'S(r,r') 



-r^G(r',r;z„)i(^.) 



+f3*(r')*Hr) 



(E3) 



Finally, the first integral in Eq. Q is transformed with 
Eq. (|16|) into an expression without K . Substituting the 
result together with Eqs. IJE2|I and (|E3p into Eq. Q, we 
arrive at Eq. (|31|1 . 



We next derive Eq. H32(l by closely following the pro- 
cedure of Ref. 1^. We first transform all the summations 
over n in Eq. Q) into contour integrals of using n{z). 
With the properties Sfl/S'i/* = and Sn/6G = 6, the 
differentiation —dfl/dT need be carried out only with 
respect to the explicit T dependence in n{z). The term 
from $ cancels that from the second term in the square 
bracket of Eq. Q). The remaining term is L of Eq. (JE1|I . 
whose explicit T dependence is transformed into the first 
term in Eq. (|E2|I . We hence obtain Eq. (|22l • 
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